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Abstract 

All of the thermodynamic information on a statistical mechanical system is encoded in the locus and density of its partition 
function zeroes. Recently, a new technique was developed which enables the extraction of the latter using finite-size data of 
the type typically garnered from a computational approach. Here that method is extended to deal with more general cases. 
Other critical points of a type which appear in many models are also studied. 
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1. Introduction 

Phase transitions are of central interest in statisti- 
cal physics and related fields. Second-order transitions 
are signaled by divergences and characterised by crit- 
ical exponents (e.g., a for the specific heat and v for 
the correlation length in the temperature driven case). 
Such non-analytic behaviour is only present in sys- 
tems of infinite extent and is therefore inaccessible to 
Monte Carlo simulations, which are restricted to a fi- 
nite number of degrees of freedom. 

Traditionally, however, finite-size scaling (FSS) 
may be used to extract thermodynamic information 
from such systems. FSS, based on the hypothesis 
that there are only two relevant length scales (namely 
the correlation length of the infinite system and its 
finite-size counterpart), typically only allows determi- 
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nation of ratios of critical exponents associated with 
thermodynamic functions, such as a/ v. An exception 
is the correlation length critical exponent v which 
can be directly extracted from logarithmic derivatives 
of magnetization moments or from the slope of the 
Binder parameter ~ (m 2k )/(m k ) 2 , k — 1,2,.... 
Alternatively one may also consider the scaling be- 
haviour of pseudocritical points. The latter, defined 
as the extrema of thermodynamic functions, approach 
the transition point as L~ x , where L denotes the lin- 
ear extent of the system and A is the so-called shift 
exponent. The exponent A coincides with 1/v in many 
models, but this is not a consequence of FSS and 
is not always true. See, e.g., [1] for a review of the 
recent literature concerning this point. A further com- 
plication that arises from the latter approach is that 
such a fit involves three parameters and is non-linear, 
so usually is quite unstable and often inaccurate. 

An increasingly popular approach is the use of FSS 
of the zeroes of the partition function. FSS of the low- 



Preprint submitted to Elsevier Science 



6 February 2008 



est zeroes in the complex temperature plane (Fisher 
zeroes) provides a direct and accurate method to ex- 
tract the exponent v, and the imaginary parts of the 
lowest zeroes (labelled by an index j) scale with lattice 
extent as Imzj ~ L~ x l v . The real part of the lowest 
partition function zero is another pseudocritical point, 
generally scaling as L~ A . 

It has long been known that a full understanding of 
the properties of the bulk system requires knowledge 
of the density of zeroes too. Until recently, determina- 
tion of the density from finite-size Monte Carlo data 
was considered difficult if not impossible. The source 
of the difficulties is that it involves reconstruction of 
a continuous density function from a discrete data set 
as the density of zeroes for a finite system essentially 
consists of a set of delta functions. 

Recent considerations have bypassed these difficul- 
ties by focusing instead on the integrated density of 
zeroes [2]. In particular, this new approach facilitates 
measurement of the strength of the transition through 
direct determination of a (as opposed to traditional 
FSS measurements of the ratio ajv). While the new 
technique proved successful, it was limited to systems 
where the zeroes fall on curves in the complex param- 
eter plane and where the zeroes are non-degenerate. 
While these two properties are common to most mod- 
els in statistical physics, they are not generic and a host 
of examples now exist where the zeroes are distributed 
across a two-dimensional region and/or occur in de- 
generate sets. Here, the new technique is extended to 
deal with such general distributions of zeroes [3]. 



2. General Distributions of Zeroes 

When the partition function, Zy, for a system 
of volume V = L d (d being the dimensionality 
of the system) can be written as a polynomial in 
an appropriate function, z, of temperature, field 
or of a coupling parameter, one has Zy(z) oc 
Y\j (z — Zj(V)) where j labels the zeroes. In the 
general case where the distribution of zeroes is two- 
dimensional, the free energy may be expressed as 
fv = JJ gv(x, y) In (z - z c - x - iy)dxdy, where 
gv is the density of zeroes and (x, y) give their loca- 
tion in the complex plane with the critical point, z c , 
as the origin. In the infinite-volume case, Stephenson 



has shown that the density near a second-order tran- 
sition point satisfies a certain homogeneous partial 
differential equation, the solution of which may be 
written as goo(x, y) = y 1 ~ a ~ m f(x/y m ), where m is 
related to the shape of the locus [4]. Integrating out 
the x— direction, and integrating up to a point r in the 
y— direction gives the cumulative density there to be 

GocMccr 2 -. (1) 

From this expression, the exponent a may be directly 
measured provided that a sensible definition for the 
cumulative density of zeroes can be applied to a finite 
system. Such a function is defined as follows. If the j th 
zero is n-fold degenerate the densities to its immediate 
left and right are given by VGy(r) = j — 1 and 
j + n — 1 respectively. The density at the j th zero, rj, 
is then defined as an average: 

G i (r,) = i(i + = -l). (2) 

Combined with (1), this allows direct determination 
of the critical exponent a. 



3. Applications 

We apply the new technique to Ising models in two 
dimensions for which the zeroes are calculable and 
which possess each of the new features we wish to 
encapsulate. In each case, the real, physical, critical 
point is characterized by a = 0. 

Brascamp-Kunz lattice with anisotropic cou- 
plings: The finite-size, standard, nearest-neighbour, 
square lattice Ising model has been solved in two 
dimensions for certain sets of boundary conditions 
including those first studied by Brascamp and Kunz 
[5]. There, periodic boundary conditions are used in 
one direction, while at the extremities in the other 
direction the spins are fixed to + on the one hand 

and the alternating sequence H 1 ... on the 

other. For a lattice of linear extent L with anisotropic 
couplings J and J', the partition function then takes 
the form of a single product (as opposed to a sum of 
four such products, which is the case when periodic 
boundary conditions in both directions are used [6]), 
greatly ameliorating the computation of its zeroes 
[1]. The zeroes are easily determined numerically 
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Fig. 1. The zeroes (top) for the anisotropic L = 40 Brascamp-Kunz 
Ising model with J'/J = 3 and their density distribution (bottom) 
near z = 1 for L = 40-140 and j = 1 (x), j = 2 (+), j = 3 
(*), j = 4 (o), j = 5 (□), j = 6 (o), j = 7 (.) and j = 8 (*). 

and are distributed across a two-dimensional region 
in the z = 2sinh(2/3) plane as shown in Fig. 1 for 
J'/J = 3 (with (3 = J/ksT). The zeroes impact 
onto the real axis at the point z = 1 and the critical 
behaviour is dominated by the zeroes close by. The 
cumulative density distribution for this set of zeroes 
is also plotted in Fig. 1. That the curve goes through 
the origin indicates the presence of a transition and 
an appropriate fit yields a = —0.016(32) compatible 
with expectations. 

Bathroom- tile lattice: Two-dimensional distribu- 
tion of zeroes may also be obtained from systems with 
isotropic couplings as demonstrated in [7] where the 
system is described in detail. In principle the full finite- 
size partition function is a sum of four product terms. 
One may construct Brascamp-Kunz type boundary 
conditions for this lattice, which would have the effect 
of projecting out one of these terms in the expression 
for the partition function. Alternatively, and more con- 
veniently, one may discretize one of the terms in the 
partition function with periodic boundary conditions 
and assume for the purposes of the analysis herein 
that the scaling behaviour of that term is generic. In- 
deed, since we are essentially interested in testing the 
scaling of the cumulative density of zeroes rather than 
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Fig. 2. The zeroes (top) for the L = 40 bathroom-tile Ising model 
and their density (bottom) near the ferromagnetic critical point 
from data with L = 40-200 and j = 1-4 (x), j = 5-12 (+), 
j = 13-16 (*), j = 17-24 (o), j = 25-32 (□), j = 33-40 (o) 
and j = 41-44 (.). 

formulating the finite lattice models themselves, this 
is sufficient for our purposes (see also [3]). 

The zeroes of such a term have varying degrees 
of degeneracy and are depicted in the complex z' = 
exp(-2/3) plane in Fig. 2, where AFM, PM, FM 
and Oi indicate the anti-ferromagnetic, paramagnetic, 
ferromagnetic and unphysical phases, respectively. 
The physical ferromagnetic critical point is given 
by z' = 0.249 038 4 . . . and the cumulative den- 
sity of zeroes nearby is also depicted in the figure. 
A fit yields a = 0.002(18), consistent with zero. 
There is also an antiferromagnetic transition point 
at z' = 4.015 445 4.... A density fit to the zeroes 
nearby yields a = 0.0006(163), again compatible 
with a = 0. 

Since the finite-size partition function is known ex- 
actly in this case, and is a convenient single product, 
it is possible to analytically extract the v and A expo- 
nents from conventional FSS of the lowest lying ze- 
roes. Indeed, one finds that v = 1 (which again gives 
a = through hyperscaling) and A = 2 in both the 
ferromagnetic and antiferromagnetic cases. 

Complex vertices: 

It has been pointed out that unphysical singular 
points (i.e., points for which there is no real (3) may 
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Fig. 3. Finite-size dependency of the zeroes close to the complex 
singular point z = i in the isotropic Brascamp-Kunz model. 
The dotted lines scale as Rez,, ~ L~ z and the solid lines as 
Rezj ~ L~ 2 . 

be considered as ordinary critical points with distinc- 
tive critical exponents (see, e.g., [7] and references 
therein). 

For the example of the two-dimensional Ising model 
with Brascamp-Kunz boundary conditions, the com- 
plex vertex points at z = i and 2i in Fig. 1 are thus 
also of interest. However, analysis shows these not to 
be of the conventional type. Fig. 3 is a plot showing 
how the first few zeroes in the simpler isotropic model 
approach such a vertex at z = i with increasing lattice 
size. 

To understand the unusual behaviour depicted, an 
analytic approach is required. A full scaling analysis of 
this model at its real transition point z = 1 is given in 
[1]. For a square lattice of linear extent L, the zeroes, 
which are labelled Zij, are given by 

Rezij — - (cos 8i + cos ) , 

Imzy- = ^Jl- (Rczij) 2 , (3) 

where 9i = (2i — l)n/L and <j>j = jn/(L + 1). Ex- 
panding the cosines for large L gives Rez^ = 1 + 
0(L~ 2 ) and Imzij = C(L _1 ), recovering A = 2 and 
v = 1, as above. For Rezy close to 1 (the physical 
critical point), both i and j have to be close to zero and 
the above expansion is legitimate. However, close to 
the unphysical point z = i, the two cosines in (3) have 
to cancel and the expansion is no longer valid. Can- 
cellations of this type remove the leading L~ 2 term 
for the approach to the vertex of the first few zeroes 
(which then scale as L~ 3 ) and do not occur at a real 
critical point. This explains the odd scaling behaviour 
in Fig. 3 and demonstrates the dangers inherent to a 



restricted traditional analysis of leading zeroes. This 
and related issues will be elaborated upon elsewhere. 

The situation close to the unphysical points z' — 
-0.601 231 8 . . . and z 1 = -1.663 251 9 ... in Fig. 2 
is more conventional and density analyses yield a = 
0.007(12) and a = -0.0095(123) respectively. The 
exponents v and A may also be extracted analytically 
here and one again finds v = 1, A = 2. It is interesting 
to note that these are yet more cases where v does not 
coincide with A [1]. 
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